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The  low-frequency  target  strength  of  shoaling  Atlantic  herring  ( Clupea  harengus )  in  the  Gulf  of 
Maine  during  Autumn  2006  spawning  season  is  estimated  from  experimental  data  acquired 
simultaneously  at  multiple  frequencies  in  the  300-1200  Hz  range  using  (1)  a  low-frequency  ocean 
acoustic  waveguide  remote  sensing  (OAWRS)  system,  (2)  areal  population  density  calibration  with 
several  conventional  fish  finding  sonar  (CFFS)  systems,  and  (3)  low-frequency  transmission  loss 
measurements.  The  OAWRS  system’s  instantaneous  imaging  diameter  of  100  km  and  regular 
updating  enabled  unaliased  monitoring  of  fish  populations  over  ecosystem  scales  including  shoals  of 
Atlantic  herring  containing  hundreds  of  millions  of  individuals,  as  confirmed  by  concurrent  trawl 
and  CFFS  sampling.  High  spatial-temporal  coregistration  was  found  between  herring  shoals  imaged 
by  OAWRS  and  concurrent  CFFS  line-transects,  which  also  provided  fish  depth  distributions.  The 
mean  scattering  cross-section  of  an  individual  shoaling  herring  is  found  to  consistently  exhibit  a 
strong,  roughly  20  dB/octave  roll-off  with  decreasing  frequency  in  the  range  of  the  OAWRS  survey 
over  all  days  of  the  roughly  2-week  experiment,  consistent  with  the  steep  roll-offs  expected  for 
sub-resonance  scattering  from  fish  with  air-filled  swimbladders. 
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I.  INTRODUCTION 

To  study  the  scattering  characteristics,  abundance  and 
diurnal  behavior  of  Atlantic  herring  ( Clupea  harengus),  the 
most  abundant  fish  species  in  and  around  Georges  Bank  dur¬ 
ing  their  Autumn  spawning  season,1  an  experiment  using 
ocean  acoustic  waveguide  remote  sensing  (OAWRS)  was 
conducted  in  the  Gulf  of  Maine  from  September  19  to  Octo¬ 
ber  6,  2006,  concentrating  on  areas  where  herring  shoals 
were  most  likely  to  form  (Fig.  1).  The  experiment,  known  as 
OAWRS  2006,  was  conducted  in  conjunction  with  the  U.S. 
National  Marine  Fisheries  Service  Annual  Atlantic  Herring 
Acoustic  Survey  of  the  Gulf  of  Maine  and  Georges  Bank. 
Fish  populations  were  instantaneously  imaged  over  a  100  km 
diameter  area  by  a  mobile  OAWRS  system”  with  minute-to- 
minute  updates  to  form  wide-area  movies  of  fish  activity 
over  many  diurnal  cycles,  demonstrating  the  capacity  of 


OAWRS  to  instantaneously  image  fish  populations  over  wide 
areas  in  complex  continental-shelf  environments  with  highly 
variable  bathymetry  and  oceanography.  Shoals  imaged  by 
OAWRS  typically  comprised  tens  to  hundreds  of  millions  of 
individuals  and  stretched  for  many  kilometers  along  the 
northern  flank  of  Georges  Bank.  Concurrent  conventional 
fish  finding  sonar  (CFFS)  surveys  showed  high  spatio- 
temporal  coregistration  with  fish  shoals  imaged  by  OAWRS 
and  provided  local  areal  population  densities,  as  well  as 
depth  and  length  distributions  of  the  fish  populations.  Con¬ 
current  trawl  sampling,  which  showed  Atlantic  herring  to  be 
the  overwhelmingly  predominant  species  comprising  the 
large  shoals,””  enabled  onsite  species  identification  and  di¬ 
rect  biological  measurements  of  parameters  such  as  fish 
length,  swimbladder  geometry,  stomach  content,  and  sexual 
development. 
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FIG.  1.  Location  of  OAWRS  2006  experiment  on  the  northern  flank  of  Georges  Bank  in  the  Gulf  of  Maine.  Plus  indicates  location  of  moored  OAWRS  source 
array  deployed  on  Oct  1-3  at  42.2089N,  67.6892W,  the  coordinate  origin  for  all  OAWRS  images  in  this  paper.  Circle  shows  typical  area  imaged  by  OAWRS, 
100  km  diameter  and  wider  than  Cape  Cod,  in  70  s.  Geographic  locations  of  trawls  deployed  by  NOAA  FRV  Delaware  11  are  overlain.  Dots  indicate  trawls 
where  herring  were  predominant  species.  In  contrast,  diamond  indicates  a  trawl  where  silver  hake  and  squids  dominated.  The  gray  dashed  box  bounds  the  area 
of  OAWRS  imaging  during  the  OAWRS  2006  experiment. 


Together  with  long-range  transmission  loss  measure¬ 
ments,  concurrent  CFFS,  and  trawl  data,  the  OAWRS  imag¬ 
ery  enabled  (1)  estimates  of  herring  target  strength  to  be 
made  at  low  frequencies  (300-1200  Hz)  from  which  physi¬ 
cal  scattering  mechanisms  may  be  inferred,  (2)  herring  spa¬ 
tial  distributions  and  abundance  to  be  estimated  over  ecosys¬ 
tem  scales,  and  (3)  regular  diurnal  patterns  in  herring 
behavior  to  be  determined  during  the  Autumn  spawning  sea¬ 
son  on  Georges  Bank.3'4 

The  mean  scattering  cross-section  of  an  individual 
shoaling  herring  is  found  to  consistently  exhibit  a  strong, 
roughly  20  dB/octave  roll-off  with  decreasing  frequency  in 
the  range  of  the  OAWRS  survey  over  many  measurement 
days,  consistent  with  the  steep  roll-off  expected  for  sub¬ 
resonance  scattering  from  fish  with  air-filled  swimbladders. 
These  findings  suggest  that  OAWRS  can  provide  valuable 
evidence  for  remote  species  classification  over  wide  areas 
since  significant  variations  in  the  frequency  dependence  of 
target  strength  are  expected  across  species  due  to  differences 
in  resonance.  This  is  because  the  dominant  source  of  acoustic 
scattering  at  low  and  mid-frequencies  (less  than  10  kHz)  is 
the  air-filled  swimbladder  for  fish  that  have  swimbladders.5 
Resonance  frequencies  depend  on  swimbladder  volume, 
shape,  ambient  pressure,  and  the  effect  of  surrounding 
tissues.6'7  For  many  fish  species  of  economic  importance  in 
the  size  ranging  from  10  to  50  cm,  resonances  are  expected 
to  range  from  several  hundreds  of  hertz  to  a  few  kilohertz  A9 
Previous  experimental  investigation  of  resonance  have  been 
limited  to  small  scale  tank  measurements  with  individual  fish 
out  of  their  natural  environment  or  highly  localized  measure  - 
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ments  made  in  situ  with  sources  of  rapid  power  roll-off  be¬ 
low  1.7  kHz.7’10’11 

The  population  that  spawns  on  the  northern  flank  of 
Georges  Bank  (Fig.  1)  is  the  largest  herring  stock  in  the  Gulf 
of  Maine,  and  has  both  ecological  and  economic 
importance.1  It  has  been  surveyed  annually  by  the  U.S.  Na¬ 
tional  Marine  Fisheries  Service  for  roughly  1  decade  during 
the  Autumn  spawning  season.1"-14  Current  estimates  of  the 
Georges  Bank  herring  stock  varies  from  500  000  to  1  X  106 
tons  based  on  acoustic  surveys  and  other  assessment 
methods,14  respectively.  The  National  Marine  Fisheries  Ser¬ 
vice  acoustic  survey  employs  highly  localized  CFFS  mea¬ 
surements  along  widely  spaced  line-transects,  roughly 
20-30  km  apart,14  trawl  sampling  at  selected  locations,  and 
takes  roughly  1  week  to  cover  the  northern  flank  of  Georges 
Bank  from  east  to  west.  As  a  result,  annual  stock  estimates 
may  be  highly  aliased  in  both  time  and  space.  One  of  the 
primary  goals  for  OAWRS  2006  is  to  provide  images  of  fish 
populations  over  the  vast  areas  they  inhabit  that  are  unaliased 
in  both  space  and  time"  so  that  more  reliable  abundance  es¬ 
timates  may  be  obtained. 

II.  MULTI-SENSOR  EXPERIMENT  DESIGN  AND 
RESOURCES 

The  OAWRS  2006  experiment  was  designed  to  coincide 
with  the  National  Marine  Fisheries  Service  annual  herring 
survey  of  Georges  Bank.  It  was  conducted  with  four  research 
vessels  (RVs)  that  employed  a  suite  of  acoustic  imaging  sen¬ 
sors,  several  oceanographic  monitoring  systems,  and  trawls 
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TABLE  I.  OAWRS  receiving  array  3-dB  angular  resolution  (i{  tf>)  at  broad¬ 
side  <f>=0  and  endfire  </>=n/2,  and  aperture  length  L  as  a  function  of  imag¬ 
ing  frequency  fc.  A  Hanning  spatial  window  is  applied  in  the  beamforming. 


fc 

(Hz) 

L 

(m) 

/3(<f>=0)l /3(</)=7t  12)  (deg) 

415 

94.5 

3.15/31.4 

735 

47.25 

3.56/33.3 

950 

47.25 

2.75/29.3 

1125 

23.625 

4.65/38.1 

for  species  identification.  The  OAWRS  vertical  source  array 
and  towed  horizontal  receiving  array  were  separately  de¬ 
ployed  from  two  medium  sized  UNOL  vessels,  RV  Endeavor 
and  RV  Oceanus  respectively,  for  bistatic  measurement  of 
echo  returns.  The  instantaneous  areal  coverage  of  the 
OAWRS  system  in  a  single  transmission  is  shown  in  Fig.  1. 
The  vertical  source  array  transmitted  a  suite  of  individual 
Tukey-windowed  linear  frequency  modulated  (LFM)  pulses 
of  1  s  duration  and  50  Hz  bandwidth  centered  at  a  suite  of 
frequencies  with  a  repetition  interval  of  150  s  for  each  center 
frequency.15  Broadband  LFM  pulses  centered  at  415  and 
735  Hz  were  transmitted  seconds  apart,  then  after  75  s  those 
centered  at  950  and  1125  Hz  were  transmitted  seconds  apart, 
and  the  process  was  repeated.  Transmissions  radiated  with 
azimuthal  symmetry  about  the  OAWRS  source  array,  for 
which  more  information  is  available  in  Ref.  16,  with  source 
level  continuously  monitored  with  two  desensitized  hydro¬ 
phones  deployed  from  RV  Endeavor.  Scattered  returns  were 
acquired  with  a  horizontal  receiving  array,  the  ONR  five- 
octave  research  array,  towed  by  RV  Oceanus  along  desig¬ 
nated  tracks.  The  multiple  nested  sub-apertures  of  the  array 
span  50-3750  Hz  frequency  range.  Returns  measured  within 
each  linear  section  of  the  array  are  processed  by  beamform¬ 
ing  and  matched  filtering  with  angular  resolution  shown  in 
Table  I.  The  receiving  array  also  contained  one  desensitized 
hydrophone  which  was  used  to  measure  transmitted  signals 
from  the  source  array  for  transmission  loss  and  source  level 
calibrations.  Two  calibrated  acoustic  targets  made  of  air- 
filled  rubber  hose,  approximately  30  m  long  and  7  cm  in 
diameter  with  known  scattering  properties,  were  vertically 
deployed  at  selected  locations  to  enable  accurate  charting  of 
scattered  returns  in  both  range  and  azimuth.  One  of  the  tar¬ 
gets  was  moored  with  lower  end  5  m  off  the  seafloor  and  the 
other  was  centered  at  140  m  in  waters  200  m  deep. 

Over  the  course  of  the  OAWRS  2006  experiment,  more 
than  3000  wide-area  images  of  the  ocean  environment  were 
acquired  by  the  OAWRS  system  for  each  of  the  four  LFM 
center  frequencies  leading  to  more  than  12  000  images  in 
total.  Similarly,  more  than  12  000  transmission  loss  measure¬ 
ments  were  made  over  the  survey  area  to  calibrate  our  trans¬ 
mission  loss  model.  The  length  of  each  RV  Oceanus  towed 
array  track  was  typically  15  km.  With  a  nominal  tow  speed 
of  2  m  s_I  for  the  receiver  ship,  a  total  of  roughly  75  images 
of  the  ocean  environment  were  generated  per  frequency 
along  each  track.  Minute-to-minute  updates  of  the  OAWRS 
imagery  made  it  possible  to  closely  monitor  herring  activity 
over  wide  areas  and  observe  patterns  of  spatial  distribution 
evolve  over  the  course  of  each  day.  The  inherent  left-right 
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FIG.  2.  (Color  online)  [(A)-(C)]  OAWRS  images  of  areal  fish  density 
zoomed-in  around  massive  herring  shoals,  with  densities  exceeding 
10  fish  m-2  in  population  centers.  Measured  during  evening  to  midnight 
hours  of  October  4,  2,  and  1,  respectively.  (A)  The  total  population  of 
herring  in  the  large  dense  shoal  is  roughly  170X  106,  and  that  in  the  diffuse 
cloud  outside  the  large  shoal  is  roughly  70  X  106.  Imaged  shoal  populations 
of  herring  are  approximately  86  X  106  and  70  X  106  respectively  for  (B)  and 
(C).  Uncertainty  in  the  abundance  estimate  is  17-20%.  Note  that  the  figures 
are  plotted  on  different  scales,  and  the  coordinate  origin  is  the  source  loca¬ 
tion  shown  in  Fig.  1. 


ambiguity  about  the  horizontal  line-array’s  axis  in  the 
OAWRS  images  were  resolved  mainly  by  varying  receiver 
ship  heading,  sometimes  only  slightly  for  several  transmis¬ 
sions  by  what  we  call  a  “Crazy  Ivan”  for  immediate  results, 
as  well  as  ship  position.  These  approaches  for  ambiguity 
resolution  are  described  in  Refs.  19-21. 

Examples  of  the  massive  fish  shoals  instantaneously  im¬ 
aged  by  OAWRS  near  Georges  Bank  are  shown  in  Fig.  2. 
The  massive  shoal  imaged  during  midnight  hours  of  October 
4  within  the  150-180  m  bathymetric  contour  in  Fig.  2(A), 
for  example,  extends  15X5  km2  and  comprises  roughly 
170  X  106  fish  distributed  about  several  population  centers. 
The  area  occupied  by  this  shoal  is  approximately  equal  to 
that  of  Manhattan  Island  in  New  York.  The  fish  population  in 
the  diffuse  cloud  region  to  the  north  is  comprised  of  over 
70  X  106  individuals. 

Concurrent  localized  imaging  of  fish  aggregations  at 
OAWRS-directed  locations  was  conducted  by  two  other  re- 
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TABLE  II.  Conventional  fish  finding  sonars,  SIMRAD  EK60  and  EK500 
specifications.  The  angular  3-dB  beamwidth  is  denoted  by  f3,  the  pulse  du¬ 
ration  by  PD,  and  repetition  rate  by  RR.  The  resolution  diameter.  Res,  is 
calculated  for  200  m  water  depth. 


Sensor 

/ 

(kHz) 

P 

(deg) 

Res 

(m) 

PD 

(ms) 

RR 

(s"1) 

EK60 

38 

120 

7 

24 

1 

1 

200 

EK500 

18 

11 

39 

2 

38 

12 

42 

1 

0.5 

120 

7 

24 

1 

search  vessels,  the  RV  Hugh  Sharp  and  the  NOAA  FRV 
Delaware  II,  using  two  downward-directed  CFFS  systems, 
the  SIMRAD  EK60  and  EK500  echosounders,  respectively. 
Both  the  EK60  and  EK500  echosounders  insonify  the  water 
column  directly  beneath  the  survey  vessel  simultaneously  at 
three  frequencies  to  provide  the  local  depth  dependence  of 
dominant  fish  layers  within  their  instantaneous  resolution 
footprints,  of  between  24-50  m  diameter,  and  estimates  of 
volumetric  and  areal  fish  population  densities.  Specifications 
of  these  two  echosounders  appear  in  Table  II.  A  Reson  7125 
Seabat  multi-beam  sonar  (400  kHz)  system  was  also  de¬ 
ployed  from  RV  Hugh  Sharp  with  an  angular  swath  of  128°. 
It  was  useful  in  providing  detailed  three-dimensional  mor¬ 
phology  of  smaller  fish  groups  located  in  the  mid-water 
column.""  A  high-speed  rope  trawl"  deployed  by  NOAA 
FRV  Delaware  II  enabled  species  identification14  at 
OAWRS-directed  locations. 

Physical  oceanography  was  monitored  by  sampling 
water-column  temperature  and  salinity  with  expendable 
bathythermographs  (XBTs)  and  conductivity-temperature- 
depth  (CTD)  sensors  at  regular  hourly  intervals  from  all  four 
research  vessels.  The  water-column  sound  speed  profile  was 
found  to  be  relatively  constant  in  space  and  time  over  the 
2006  OAWRS  survey,  as  shown  by  the  compilation  of  over 
roughly  200  samples  taken  during  the  experiment  in  Fig.  3. 
No  mesoscale  oceanographic  features  such  as  eddies  were 
found  or  expected.  The  small  fluctuations  about  the  mean 
profile  are  due  to  mild  internal  wave  activity  that  causes 
well-understood  short-term  Gaussian  field  fluctuations  in 
acoustic  transmission  that  have  an  intensity  standard  devia¬ 
tion  that  can  be  reduced  to  a  small  fraction  of  the  mean  by 
stationary  averaging. 24-26  An  instrumented  tow  cable  was 
also  deployed  from  the  RV  Hugh  Sharp  to  provide  continu¬ 
ous  measurement  of  temperature.  This  oceanographic  infor¬ 
mation  was  used  to  carefully  update  horizontal  locations 
and  depths  of  the  OAWRS  source  (typically  centered  at 
60-70  m)  and  receiving  arrays  (centered  at  105  m)  (Ref.  15) 
to  optimize  OAWRS  imaging  of  fish  groups.  Decisions  were 
often  based  on  the  outputs  of  the  range-dependent  acoustic 
model  (RAM),  based  on  the  parabolic  equation,  for  multi¬ 
modal  waveguide  transmission  loss  in  the  range-dependent 
Georges  Bank  environment  with  the  hourly  sound  speed  pro- 
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FIG.  3.  Profiles  of  water-column  sound  speed  from  XBT  and  CTD  measure¬ 
ments  made  from  all  four  research  vessels  on  the  Northern  Flank  of  Georges 
Bank  and  Georges  Basin  during  OAWRS  2006. 


III.  DATA  PROCESSING  AND  ANALYSIS 

A.  Generating  instantaneous  wide-area  OAWRS 
images  of  the  ocean  environment 

Wide-area  images  of  instantaneous  scattered  intensity 
spanning  100  km  in  diameter  were  generated  in  near  real¬ 
time  for  every  broadband  transmission  centered  at  each  of 
the  four  frequencies,  fc= 415,  735,  950,  and  1125  Hz.  For 
each  transmission,  the  pressure  data  on  the  receiving  array 
were  first  beamformed  to  determine  the  azimuth  of  the  arriv¬ 
als,  then  matched  filtered  with  the  source  signal,  and  charted 
in  range  using  two-way  travel  time.19-"1  Each  image  was 
then  mapped  onto  geographic  space  using  the  GPS  latitude 
and  longitude  information  of  the  source  and  receiving  array. 
A  nominal  sound  speed  of  1475  m/s  that  minimizes  charting 
errors  was  used  to  convert  the  travel-time  of  the  signal  to 
range." 1,28  The  range  resolution  Ap  of  the  OAWRS  system  is 
approximately  15  m  after  matched- filtering,  and  the  azi¬ 
muthal  resolution  /3{(f>,fc )  associated  with  each  frequency 
band  both  at  broadside  and  endfire  is  tabulated  in  Table  I.  A 
Hanning  spatial  window  was  applied  in  the  beamforming  to 
reduce  sidelobe  levels  by  more  than  30  dB  from  the  main 
lobe.  A  detailed  explanation  of  the  image  formation  process 
is  provided  in  Refs.  2,  20,  and  21. 

A  standard  procedure  of  averaging  three  consecutive  in¬ 
stantaneous  OAWRS  images  and  two  adjacent  range  cells  is 
used  for  all  OAWRS  images  presented  here.  This  leads  to  an 

experimentally  determined  standard  deviation  in  log- 

"2 

intensity  of  roughly  1.5  dB,  consistent  both  with  theory  and 
previous  experiments.  ’  ’  ’  ’  This  standard  deviation  is 
negligible  compared  to  the  dynamic  ranges  of  features  in  the 
OAWRS  images  and  the  variations  in  herring  target  strength 
measured  across  frequency  in  the  OAWRS  range. 
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B.  Estimating  areaf  fish  population  density  from 
instantaneous  OAWRS  imagery 


Here  we  describe  how  areal  fish  population  density  over 
wide  areas  may  be  estimated  from  OAWRS  intensity  images 
of  the  ocean  environment.  At  typical  OAWRS  operating  fre¬ 
quencies  from  hundreds  of  hertz  to  a  few  kilohertz,  most  fish 
are  acoustically  compact  scatterers,  with  swimbladder  sizes 
that  are  much  smaller  than  the  wavelength.  The  sonar  equa¬ 
tion  approach  is  then  valid  for  analyzing  scattering  from  fish 
since  the  scattered  field  from  each  individual  is  omni¬ 
directional,  making  propagation  and  scattering  factorable 
even  in  a  waveguide."9  The  expected  scattered  intensity  from 
fish  aggregations  after  matched-filtering  is  dominated  by  the 
incoherent  intensity  or  variance  of  the  scattered  field  and 
multiple  scattering  effects  are  negligible  for  the  densities 
found  here,  as  shown  in  Refs.  30  and  31.  As  a  result,  given  a 
source  at  r0  transmitting  a  broadband  signal  with  bandwidth 
centered  at  fc  and  a  receiver  at  r,  the  expected  scattered 
intensity,  ( Is(pm,fc )),  within  the  OAWRS  resolution  footprint 
of  area  A(pm\Ap,fc)  centered  at  horizontal  location  pm  can 
be  expressed  as 

M(pm) 

(UPmJc))=  2  (\Q(fc)\2) 

1=1 

X  (4  tt)4<  I  G(r,|r0,/C)  G(r \rhfc)  |: 2) 

x(fe^),  0, 


where  M(pm )  is  the  number  of  fish  within  the  resolution  cell, 
r,  is  the  location  of  the  z'th  fish,  |<2(/c)|2  is  the  source  inten¬ 
sity,  G(r,|r0,/C)  and  G{r\rhfc)  are  the  waveguide  Green’s 
functions  from  the  source  to  each  scatterer  and  from  each 
scatterer  to  the  receiver,  respectively,  S(rhfc)  is  the  fish  scat¬ 
ter  function,  and  k  is  the  wavenumber. 

The  expected  intensity  in  a  fluctuating  waveguide  from 
uniformly  distributed  targets  within  the  resolution  footprint 
can  be  approximated  as 


M(p,„) 

( UPmJc ))  “  (\Q(fc)\2)y(Pnvfc)  S 


i=  1 


1s(g-'/c)12\ 

e  / 


(2) 


where  y(pm ,/,.)  =  (|  (4  tt) 2C(r„,  |  r0  ,/c)  G(r  \  r,  ,/c)  | 2)  for  suffi¬ 
ciently  narrow  depth  layers  H  and  areal  footprints  over 
which  y(pm,fc)  becomes  effectively  constant,  as  shown  for 
the  OAWRS  2006  fish  shoal  imaging  in  Ref.  4.  The  last 
factor  of  Eq.  (2)  can  be  written  as 


M|”7lA(r,/c)|2\ 


|S(r,-,/c)|2 

k 2 


P(r  i)dr] 


=  M{pm)cr{pm,fc),  (3) 

where  P(r;)  is  the  probability  density  of  finding  the  z'th 
fish  at  location  rh  and  P{r,)=  1 IA(pm\Ap,fc)H  for  uniformly 
distributed  fish  shoals,  cf(pm,/c)  is  the  average  scattering 
cross-section  of  an  individual  fish  over  the  OAWRS  reso¬ 
lution  footprint  and  the  depth  layer,  nAOAWRS(pm)  = 
M(pm)l A(pm\Ap,fc)  is  the  mean  areal  fish  population  den¬ 
sity  within  the  resolution  footprint,  and  A(pnl\Ap,fc)~ 


pmAp/3{<J),fc)  is  the  range  and  azimuth-dependent  spatial 

21 

resolution  of  the  OAWRS  imaging  system. 

Inserting  Eq.  (3)  into  Eq.  (2)  and  taking  10  log10  of  both 
sides,  we  obtain  the  scattered  intensity  level  in  decibels; 

L(p,mfc)  ~  SL(fc)  +  TTL(pm,fc)  +  SSoawrs (Pm>/c) 

+  10  log10(A(pm|Ap,/c)) ,  (4) 

where  T(pm,/C)  =  101og10(/S(pm,/C)),  TTL{pm,fc)  = 
10  log loyiPnnfc)  describes  the  expected  second  moment  of 
depth  averaged  propagation  to  and  from  the  fish  layer  aver¬ 
aged  over  the  resolution  footprint  of  the  OAWRS  system, 
SL(fc)  =  1 0  log  10<  |  Q{fc)  | 2)  is  the  spectral  source  level,  and 
SS0AWRs(Pm,/c)  is  the  scattering  strength. 

From  Eqs.  (2)-(4),  OAWRS  scattering  strength  can  be 
expressed  as 

S  S  OAWRS  (Pm  >fc)  =  TS0AWRs(/c) 

+  10  logio(«A,OAWRs(Pm))>  (5) 

where  TSOAwrs(//)  =  10  \og  wa(pm,fc)  in  units  of  dB  re  1  nr 
is  the  target  strength  corresponding  to  the  average  scattering 
cross-section  of  an  individual  fish  over  the  OAWRS  reso¬ 
lution  footprint  and  depth  layer  within  the  bandwidth  cen¬ 
tered  at  fc. 

The  terms  in  Eq.  (4)  are  evaluated  separately  for  each  of 
the  four  OAWRS  LFM  waveforms  with  different  center  fre¬ 
quencies  fc.  A  calibrated  stochastic  transmission  loss  model 
based  on  the  parabolic  equation"  for  a  range-dependent  fluc¬ 
tuating  ocean  waveguide  is  used  to  estimate  the  random 
Green’s  functions  and  determine  TTL(pm,fc)  following  the 
approach  described  in  the  Appendix  E  and  Ref.  16.  Expected 
source  level  is  estimated  from  one-way  propagated  signals 
received  by  a  desensitized  hydrophone  on  the  moving  re¬ 
ceiver  array  using  the  approach  of  Ref.  16.  The  two  moni¬ 
toring  hydrophones  on  the  source  ship  were  used  to  verify 
the  source  level  estimates.  Our  analysis  indicates  the  source 
transmitted  a  stable  output  over  the  course  of  each  day. 

The  application  of  Eq.  (4)  to  estimate  scattering  strength 
from  OAWRS  imagery  is  illustrated  in  Figs.  2  and  3  of  Ref. 
32  and  in  Ref.  33.  Scattering  strength  is  a  useful  parameter 
for  characterizing  submerged  objects,  both  distant  and 
nearby,  because  it  is  independent  of  the  spatially  varying 
transmission  loss  and  areal  resolution  footprint  of  the  imag¬ 
ing  system.  Once  the  target  strength  expected  of  an  indi¬ 
vidual  fish  is  known,  an  areal  fish  population  density  image 
can  be  obtained  from  a  scattering  strength  image" "  using  Eq. 
(5).  The  target  strength  corresponding  to  the  average  scatter¬ 
ing  cross-section  of  an  individual  fish  at  OAWRS  operating 
frequencies  is  estimated  by  matching  between  OAWRS  and 
CFFS  areal  fish  population  density  measurements  where  si¬ 
multaneous  sampling  through  stationary  fish  populations  is 
available. 


C.  Estimates  of  areal  fish  population  density  from 
CFFS 

The  CFFS  measurements  at  38  kHz  are  used  to  provide 
local  estimates  of  areal  fish  population  density.14'"4"0  The  7° 
3-dB  beamwidth  yields  an  instantaneous  circular  survey  area 
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TABLE  III.  Physical  parameters  of  modeled  fish  species  and  their  measured 
target  strength  at  38  kHz  with  a  CFFS. 


Species 

Atlantic  herring 

Acadian  redfish 

Silver  hake 

Lfl  (cm) 

19-30 

15-39 

2-35 

Ltl  (cm) 

20-34 

16-41 

2-35 

Depth  (m) 

120-190 

120-190 

10-75 

<r5CFFS>  (dB) 

-39. 7a 

-38.9b 

N/A 

^CFFS  (dB) 

1.3C 

2.4 

N/A 

p/  (kg  m-3) 

1071 

1080 

1050 

f  (Pa  s) 

50 

50 

20 

f 

€ 

5-10 

6 

8 

Knb 

0.05 

0.05 

0.03 

h 

Kz 

x1 

0.05 

0.03 

JMean  target  strength  of  herring  calculated  using  Eq.  (7). 
bMean  target  strength  of  redfish  calculated  using  equation  in  Ref.  71. 
“Standard  deviation  of  derived  herring  target  strength  at  38  kHz  incorporat¬ 
ing  fish  length  and  depth  distribution  from  CFFS  and  trawl  surveys. 
dFish  flesh  density. 

“Viscosity  of  fish  flesh. 

'major-to-minor-axis  ratio  of  fish  swimbladder. 

gFish  swimbladder  volume  to  fish  body  volume  ratio  at  neutral  buoyancy 
depth. 

hFish  swimbladder  volume  to  fish  body  volume  ratio  at  depth, 
x  is  a  linear  function  of  ambient  pressure  at  depth  given  by  x 
=  Knb(Pnbl  Px),  where  Pnb  is  the  ambient  pressure  at  neutral  buoyancy  depth 
Znb,  and  Pz  is  the  ambient  pressure  at  any  depth  z. 

of  24  m  diameter  directly  under  the  survey  vessel  at  200  m 
water  depth.  Volumetric  scattering  from  all  targets  within  the 
conical  beam  were  measured.  The  localized  areal  fish  popu¬ 
lation  density  in  fish  m“2,  denoted  by  n4  CHhS,  can  be  esti¬ 
mated  using 

477  fZl 

ma,cffs  =  _  I  svdz,  (6) 

°bsd  zj 

where  sv  is  the  volume  backscattering  coefficient36  in  m-1,  z,\ 
and  7.2  delimit  the  depth  bounds  for  fish  aggregations,  and 
cfbs=47rlC/TScFFS^/10  is  the  mean  backscattering  cross-section 
of  an  individual  at  38  kHz  in  units  of  m2,  where  (TSChbs)  is 
the  corresponding  mean  target  strength  at  ultrasonic  fre¬ 
quency. 

The  expected  target  strength  for  an  individual  fish  at 
38  kHz  varies  with  species,  depth,  and  total  fish  length. 
Here,  the  expected  TScffs  in  dB  of  an  individual  herring  of 
total  length  Ln  in  centimeters  at  depth  z  in  m  is  obtained 
from  Eq.  (5)  of  Ref.  37, 

TScffs  =  20  log10  Lxl  -  2.3  log(l  +  zj  10)  -  65.4,  (7) 

and  then  converted  to  crbs.  The  mean  backscattering  cross- 
section  (rbs  is  obtained  as  the  weighted  average  over  the  total 
length  and  depth  distribution  of  the  fish  aggregations.  From 
trawl  surveys  of  the  imaged  fish  populations  in  OAWRS, 
herring  was  the  overwhelmingly  dominant  species  compris¬ 
ing  the  large  shoals,  which  had  small  fractions  of  redfish  and 
silver  hake.  Estimates  of  the  mean  TScffs  for  individual  her¬ 
ring  and  redfish  based  on  our  trawl  measurements  (Appendix 
A)  of  the  length  distribution  are  provided  in  Table  III.  The 
expected  target  strength  of  herring  and  redfish  over  similar 
depth  extent  at  38  kHz  are  close,  varying  at  most  by  1  dB, 
albeit  their  different  length  distributions.  In  contrast,  their 
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low-frequency  target  strength  near  resonance  varies  signifi¬ 
cantly,  as  discussed  in  Sec.  IV  C.  These  modeled  TScffs  val¬ 
ues  are  in  good  agreement  with  those  obtained  by  experi¬ 
mentally  analyzing  the  CFFS  backscattered  field  from 
individual  fish  distinguishable  in  the  periphery  of  various 
aggregations  consistent  with  0.1  dB  mean  squared  errors  re¬ 
ported  in  Ref.  37. 

Figures  4(D)  and  4(E)  illustrate  the  application  of  Eq. 
(6)  to  estimate  areal  fish  density  for  herring  aggregations  in 
the  120-180  m  water  depth  range. 

D.  Estimating  low-frequency  target  strength  by 
matching  OAWRS  and  CFFS  population  densities 

Here  we  describe  our  procedure  for  estimating  the  low- 
frequency  target  strength  corresponding  to  the  average  scat¬ 
tering  cross-section  of  an  individual  shoaling  herring  over 
the  resolution  footprint  of  the  OAWRS  system  by  correlating 
OAWRS  data  with  simultaneous  measurements  made  along 
CFFS  transects.  The  target  strength  of  herring  at  950  and 
1125  Hz  is  found  to  be  significantly  higher  than  at  415  and 
735  Hz,  making  much  lower  herring  densities  observable  at 
these  higher  frequencies.  At  the  lower  frequencies  of  415  and 
735  Hz,  the  herring  target  strength  is  weaker  causing  the 
scattered  returns  to  be  background  saturated  at  moderate  fish 
densities.  Due  to  the  receiving  array’s  sub-aperture  design, 
OAWRS  images  at  950  Hz  have  the  best  cross-range  reso¬ 
lution,  making  this  an  optimal  frequency  for  wide-area  sens¬ 
ing.  An  alternative  approach  for  target  strength  estimation, 
based  on  differencing  pairs  of  OAWRS  wide-area  scattering 
strength  images  at  two  distinct  frequencies,  is  applied  in  Sec. 
Ill  E  to  determine  target  strength  at  the  lower  frequencies. 
The  target  strength  estimates  are  summarized  in  Table  IV. 

Calibrated  acoustic  targets  were  deployed  on  October 
2-3  enabling  independent  and  precise  geographic  charting  of 
OAWRS  images.  By  making  small  adjustments  to  the  chart¬ 
ing  speed  and  array  orientation,  scattered  returns  from  cali¬ 
brated  targets  were  accurately  charted  to  the  correct  range- 
azimuth  resolution  cell  relative  to  the  source  and  receiver. 
This  ensures  that  scattered  returns  from  all  other  targets,  in¬ 
cluding  the  fish  aggregations,  have  been  accurately  charted 
as  well.  In  this  section  we  focus  on  data  acquired  on  October 
2-3  when  calibrated  target  data  were  available  and  present 
target  strength  estimates  for  other  days  in  the  Appendix  B. 

Close  to  midnight  on  October  2,  both  OAWRS  and 
CFFS  systems  simultaneously  co-registered  a  massive  her¬ 
ring  shoal  between  the  150  and  180  m  isobaths  on  the  north¬ 
ern  flank  of  Georges  Bank,  as  shown  in  Figs.  4(A)-4(C).  The 
observations  were  made  continuously  over  a  90-min  period 
between  23:30  Eastern  Daylight  Time  (EDT)  on  October  2 
and  01:00  EDT  on  October  3.  Measurements  from  the  two 
systems  are  highly  correlated  during  the  course  of  the  obser¬ 
vations  because  of  the  statistical  stationarity  of  the  fish  popu¬ 
lations  even  though  their  resolution  footprints  are  signifi¬ 
cantly  different.  The  OAWRS  system  monitored  and  sampled 
the  temporal  and  spatial  evolution  of  the  shoal’s  horizontal 
morphology  at  intervals  of  75  s  without  aliasing.  Concur¬ 
rently,  the  CFFS  system  crossed  the  same  shoal  twice  along 
a  U-shaped  transect  with  two  parallel  transects  1.5  km  apart. 
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FIG.  4.  (Color  online)  Herring  target  strength  at  950  and  1125  Hz  estimated  by  matching  areal  fish  density  in  OAWRS  and  CFFS  data  acquired  during 
midnight  hours  of  October  2.  [(A)-(C)]  A  sequence  of  instantaneous  OAWRS  scattering  strength  images  zoomed  into  the  region  containing  a  massive  herring 
shoal  with  overlain  CFFS  line-transect  (solid  line)  made  at  nominal  tow-speed  of  2.5  m  s-1.  (D)  CFFS  time-depth  echogram  provides  local  depth  distributions 
of  fish  aggregations.  Dashed  lines  at  23:30  EDT  and  01 :00  EDT  correspond  to  transect  start  and  end  points  a  and  11.  respectively.  (E)  The  areal  fish  population 
densities  inferred  from  CFFS  measurements  following  Eq.  (6)  are  plotted  as  a  function  of  time  in  black,  and  the  corresponding  areal  fish  population  densities 
in  dB,  10  log10(n/i  cffs)>  are  plotted  in  gray.  (F)  The  OAWRS  scattering  strength  measurements  and  (G)  instantaneous  target  strength  estimates  along  CFFS 
line-transects  at  950  and  1125  Hz.  Target  strength  estimates  near  the  edge  of  shoals  are  not  accurate  because  of  nonstationarity.  (H)  Population  of  herring 
within  the  area  shown  in  (A)-(C)  determined  with  various  OAWRS  fish  density  nA  thresholds.  Solid  line  gives  population  above  the  threshold  and  dotted  line 
gives  population  below  the  threshold. 


The  depth  distribution  of  the  fish  population,  within  roughly 
40-60  m  of  the  seafloor,  is  relatively  consistent  across  the 
two  CFFS  transects,  as  shown  in  Fig.  4(D). 

To  accurately  estimate  low-frequency  target  strength,  we 
confine  our  present  analysis  to  contiguous  space-time  seg¬ 
ments  that  consistently  register  significant,  stationary  scatter¬ 
ing  from  fish  aggregations  in  both  OAWRS  and  CFFS  sys- 
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terns.  We  derive  threshold  values  for  CFFS  population 
density  and  OAWRS  scattering  strength.  The  segmented  data 
above  these  thresholds  are  used  for  target  strength  estima¬ 
tion.  The  CFFS  threshold  is  set  at  0.2  fish  nr2,  as  shown 
in  Fig.  4(E).  For  OAWRS,  two  square  areas  of  dimension 
6.2X3  and  1.57X5.58  km2  that  continuously  register  sig¬ 
nificant  fish  scattering  and  diffuse  background  reverberation, 
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TABLE  IV.  Mean  low-frequency  target  strength  estimates.  The  TScor  esti¬ 
mates  are  obtained  by  correlating  OAWRS  with  CFFS  data  along  CFFS 
transect.  This  approach  is  only  applied  to  OAWRS  data  at  950  and  1 125  Hz. 
For  the  other  frequencies,  the  TSsc  estimates  are  obtained  by  the  approach  of 
differencing  OAWRS  images.  The  Diff  is  the  expected  target  strength  dif¬ 
ference  between  the  given  frequency  and  950  Hz. 


Date 

fc 

(Hz) 

TS 

1  °cor 

(dB  re  1  m2) 

Diff 

(dB) 

TSsc 

(dB  re  1  m2) 
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-44.4 
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1125 

N/A 

=  10 

=-34.4 

2 

“Base  target  strength  at  950  Hz  used  in  the  differencing  approach. 


respectively,  throughout  the  course  of  observation  are  first 
examined.  The  histogram  of  scattering  strength  values  within 
these  areas,  averaged  over  multiple  OAWRS  images,  are 
plotted  in  Fig.  5(D).  The  histograms  are  approximately 
Gaussian.  The  OAWRS  threshold  is  then  set  at  -50  dB  for 
950  Hz  to  distinguish  fish  scattering  from  the  background. 
This  threshold  is  roughly  2  standard  deviations  below  the 
mean  for  the  fish  histogram  and  roughly  2  standard  devia¬ 
tions  above  the  background  mean. 

Employing  Eq.  (5),  and  assuming  local  stationarity  of 
fish  population,  we  set  the  areal  fish  density  within  the 
OAWRS  resolution  footprint  to  that  simultaneously  sampled 
by  CFFS  transect  through  the  OAWRS  footprint,  ffi.oAWRS 
~nA cffs-  The  resulting  target  strength  estimates  for  fish  in 
these  contiguous  shoaling  regions  at  950  Hz  are  shown  in 
Fig.  4(G).  The  differences  in  target  strength  estimates  along 
the  transect  are  due  to  the  fact  that  the  OAWRS  and  CFFS 
systems  have  different  resolution  footprint  sizes,  and  so  the 
true  mean  fish  areal  densities  within  the  OAWRS  resolution 
cell  may  be  overestimated  or  underestimated  by  the  CFFS 
system  given  nonstationary  spatial  distributions,  as  occurred 
at  shoal  boundaries.  The  combination  of  measurements  from 
many  space-time  locations  from  both  systems  should  yield 
mean  target  strength  estimates  with  small  variance  by  virtue 
of  the  law  of  large  numbers  as  discussed  in  Appendix  C. 
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Similar  statistical  analyses  have  been  conducted  for 
OAWRS  data  at  1125  Hz,  with  estimated  target  strength  ap¬ 
pearing  in  Fig.  4(G).  This  approach  is  also  applied  to  infer 
herring  target  strength  at  950  and  1125  Hz  using  OAWRS 
and  CFFS  data  on  October  3,  where  two  contiguous  shoal 
segments  are  imaged.  Estimated  target  strengths  for  these 
two  segments  are  provided  in  Table  IV.  The  target  strength 
estimates  at  950  Hz  for  the  three  data  sets  are  consistent, 
with  a  standard  deviation  of  roughly  1  dB. 

The  approach  of  this  section  is  not  used  to  estimate  tar¬ 
get  strength  at  415  and  735  Hz  because  the  herring  are  much 
weaker  scatterers  at  these  frequencies,  as  seen  in  Fig.  5(A) 
where  only  the  densest  shoal  population  centers  stand  above 
background  scattering  levels.  An  alternative  approach  to  es¬ 
timate  herring  target  strength  at  these  lower  frequencies  is 
developed  and  applied  in  Sec.  Ill  E. 

E.  Frequency  dependence  of  target  strength 
estimated  by  differencing  OAWRS  scattering  strength 
images  over  wide  areas 

Here,  we  develop  an  alternative  approach  to  estimate 
target  strength  expected  of  an  individual  shoaling  herring  by 
differencing  pairs  of  OAWRS  scattering  strength  images  ac¬ 
quired  at  two  distinct  frequencies  over  the  entire  area  of  the 
shoal.  We  apply  this  to  data  at  415  and  735  Hz.  The  ap¬ 
proach  is  illustrated  by  the  conceptual  diagram  shown  in  Fig. 
6.  From  Eq.  (5),  we  observe  that  scattering  strength  in  areas 
containing  fish  increases  logarithmically  with  areal  density 
nA.  Here,  j\  represents  a  low  frequency,  such  as  415  Hz, 
where  the  target  strength  for  fish  is  lower  than  at  another 
frequency  f2,  such  as  950  Hz.  The  background  scattering 
strength  from  sources  other  than  fish  is  expected  to  be  statis¬ 
tically  stationary  and  can  be  identified  by  its  mean  level 
which  is  frequency  dependent.  The  total  scattering  strength 
measured  at  any  given  pixel  is  a  sum  of  the  contribution 
from  fish  and  other  background  effects.  The  difference  be¬ 
tween  the  total  scattering  strength  across  various  pixels  at  the 
two  frequencies  then  follows  the  trend  illustrated  in  Fig. 
6(B),  where  at  very  low  fish  densities,  the  scattering  strength 
is  dominated  by  the  background  reverberation,  and  at  very 
high  fish  densities  by  fish  scattering.  The  difference  in  the 
scattering  strength  at  low  fish  densities  therefore  provides  a 
measure  of  the  difference  in  background  reverberation.  The 
difference  in  the  scattering  strength  at  high  fish  densities  is 
equal  to  the  target  strength  difference  for  fish  at  these  fre¬ 
quencies.  If  the  target  strength  at  one  of  the  frequencies  is 
known  accurately,  then  the  target  strength  at  the  other  fre¬ 
quency  can  be  obtained. 

This  approach  is  implemented  for  pairs  of  OAWRS  im¬ 
ages  using  950  Hz  as  the  base  frequency.  The  difference  in 
scattering  strength  is  calculated  and  plotted  for  OAWRS  data 
acquired  between  22:00  and  22:45  EDT  on  October  3  in 
Figs.  7(A)-7(C)  for  various  frequency  pairs.  We  observe  the 
scattering  strength  difference  in  the  background  is  roughly 
1  dB  between  1125  and  950  Hz,  but  the  fish  target  strength 
difference  is  larger,  roughly  7.5  dB.  Between  950  and 
415  Hz,  the  background  scattering  strength  difference  is 
roughly  1.5  dB,  but  the  fish  target  strength  difference  is  more 
than  17  dB.  Between  950  and  735  Hz,  no  conclusion  can  be 
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FIG.  5.  (Color  online)  The  intensity  of  scattered  returns  from  shoals  is  highly  frequency-dependent.  The  histograms  illustrate  that  it  is  easier  to  detect  shoals 
over  background  regions  at  higher  frequencies.  Simultaneous  trawls  show  shoals  are  overwhelmingly  comprised  of  herring  while  background  regions  yield 
negligible  herring  (Table  IV,  Fig.  12).  [(A),  (C),  and  (E)]  OAWRS  images  of  herring  shoal  acquired  simultaneously  at  three  distinct  frequency  bands  centered 
at  415,  950,  and  1125  Hz  at  00:41:15  EDT  on  October  2.  The  colorscale  used  in  (A),  (C),  and  (E)  is  the  same  as  in  Figs.  4(A)-4(C).  [(B),  (D),  and  (F)] 
Histograms  of  scattering  strength  values  at  locations  within  the  shoal  (areas  inside  the  dashed  box)  and  in  a  background  region  (areas  inside  the  solid  box) 
plotted  for  comparison.  The  735  Hz  data  are  ambient  noise  limited  in  background  areas  due  to  weak  source  level  and  is  not  shown. 


drawn  about  background  scattering  because  the  735  Hz  data 
were  dominated  by  ambient  and  nearby  shipping  noise  since 
the  source  level  for  this  frequency  was  lower.  The  fish  target 
strength  difference  between  950  and  735  Hz  is  roughly 
7.5  dB.  These  results  are  tabulated  in  Table  IV. 

As  the  OAWRS  imaging  frequency  increases  from 
415  to  1125  Hz,  the  target  strengths  of  both  fish  populations 
and  background  levels  also  increase.  From  the  histograms  of 
Fig.  5,  the  scattering  strength  increase  with  frequency  is 
greater  for  the  fish  shoals  than  background  levels,  making  it 
easier  to  detect  fish  aggregations  at  the  higher  frequency.  At 
415  Hz,  only  densely  populated  fish  regions  are  distinguish¬ 
able  from  the  background.  Fish  densities  at  shoal  peripheries 
are  typically  too  low  to  be  detectable. 

The  background  levels  in  Figs.  7(A)-7(C)  can  be  used  to 
derive  the  minimum  detectable  fish  densities  in  the  OAWRS 
system  at  various  frequencies.  From  Fig.  5(D),  if  we  require 
that  fish  returns  stand  at  least  1  standard  deviation  above  the 
background  mean  to  be  detectable,  then  scattering  strength 
levels  that  are  above  roughly  -52  dB  at  950  Hz  would  be 
detectable.  This  corresponds  to  a  minimum  detectable  areal 
density  of  roughly  0.1  fish  in  2.  The  minimum  detectable 
densities  at  the  other  frequencies  are  tabulated  in  Table  V. 
These  are  based  on  scaling  the  fish  densities  up  or  down 
depending  on  the  background  mean  scattering  strength  level 
at  the  other  frequencies  (Fig.  7)  and  also  accounting  for  the 
target  strength  differences.  These  results  are  consistent  with 
those  obtained  from  analyzing  the  histograms  in  Fig.  5. 
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IV.  RESULTS  AND  DISCUSSION 
A.  Measured  abundance 

Herring  areal  population  densities  and  abundances  are 
estimated  by  subtracting  the  estimated  target  strength  ex¬ 
pected  of  an  individual  herring  from  OAWRS  scattering 
strength  images  as  explained  in  Sec.  Ill  B.  Areal  fish  density 
is  first  calculated  by  applying  Eq.  (5)  at  each  pixel  in  an 
OAWRS  image.  Integrating  over  the  area  of  a  shoal  then 
provides  an  estimate  of  its  abundance.  We  illustrate  abun¬ 
dance  estimation  with  OAWRS  images  generated  at  950  Hz 
since  these  have  the  best  spatial  resolution. 

Figures  4(A)-4(C)  show  the  areal  fish  density  for  a  se¬ 
quence  of  instantaneous  OAWRS  images  close  to  the  mid¬ 
night  hours  of  October  2  where  average  density  within  the 
shoal  often  exceeds  10  fish  m-2.  Figure  4(H)  shows  the 
population  over  time  within  the  area  shown  in  Figs. 
4(A)-4(C)  for  consecutive  OAWRS  images  from  23:30  EDT 
on  October  2  to  01:00  EDT  on  October  3.  When  no  threshold 
is  applied,  we  simply  integrate  the  densities  throughout  the 
area.  The  total  population  of  roughly  40  X  106  fish  includes 
fish  in  the  shoaling  region  as  well  as  diffuse  fish  clouds  out¬ 
side  the  shoal.  Contributions  from  outside  the  fish  shoal, 
which  may  include  background  from  the  seafloor,  are  esti¬ 
mated  to  account  for  less  than  10%  of  the  total  population  in 
the  area  shown. 

To  exclude  background  reverberation,  a  density  thresh¬ 
old  is  selected  to  segment  shoaling  regions  where  fish  scat- 

Gong  et  a!:.  Target  strength  and  abundance  of  shoaling  herring 


Downloaded  21  Dec  201 1  to  18.38.0.166.  Redistribution  subject  to  ASA  license  or  copyright;  see  http://asadl.org/journals/doc/ASALIB-home/info/terms.jsp 


TABLE  V.  Minimum  detectable  fish  density  (M.D.D)  in  OAWRS  imagery. 


FIG.  6.  (A)  Schematic  of  scattering  strength  levels  at  two  distinct  frequen¬ 
cies  fi  and  f2  as  a  function  of  local  fish  density  and  (B)  their  difference  at  a 
given  OAWRS  pixel.  The  scattering  strength  difference  equals  the  mean 
background  level  difference  for  low  fish  densities,  while  at  high  fish  densi¬ 
ties,  the  scattering  strength  differences  equals  the  fish  target  strength  differ¬ 
ence. 


tering  is  dominant.  Figure  4(H)  shows  the  total  population  of 
fish  above  and  below  various  density  thresholds,  0.1,  0.2, 
and  0.5  fish  m”2.  The  optimal  density  threshold  for  segment¬ 
ing  the  shoaling  region  is  roughly  0.2  fish  m  2  since  below 
this  threshold  populations  stay  fairly  constant  as  expected  for 
background  levels.  This  threshold  also  corresponds  to  a  scat¬ 
tering  strength  of  roughly  -50  dB  which  is  roughly  the  value 
where  the  scattering  strength  histograms  for  background  and 
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shoaling  regions  intersect  in  Fig.  5.  The  population  of  fish  in 
the  shoaling  region  varies  over  time  between  30  X  106  and 
40  X  106. 

The  massive  shoal  shown  in  Fig.  2(A)  is  comprised  of 
over  240  X  106  fish,  with  roughly  170  X  106  in  the  large  con¬ 
solidated  shoal  and  70  X  106  in  the  diffuse  fish  aggregation 
region.  The  population  of  fish  in  the  shoals  of  Figs.  2(B)  and 
2(C)  are  roughly  86  X  106  and  70  X  106,  respectively.  As  dis¬ 
cussed  in  the  Appendix  D,  the  uncertainty  in  the  abundance 
estimate  is  roughly  17-20%. 

B.  Measured  low-frequency  target  strength 

Following  the  approach  of  Secs.  Ill  D  and  III  E,  the  tar¬ 
get  strength  corresponding  to  the  mean  scattering  cross- 
section  of  an  individual  shoaling  herring  is  estimated  as  a 
function  of  frequency  in  the  300-1200  Hz  range  of  the 
OAWRS  system  from  measured  scattered  returns  and  mea¬ 
sured  and  modeled  transmission  loss  as  shown  in  Fig.  10  and 
Table  IV.  The  target  strength  data  show  a  consistent  depen¬ 
dence  in  both  level  and  roll-off  with  decreasing  frequency  of 
roughly  20  dB/octave  for  all  measurements,  which  spanned 
five  shoals  on  4  days.  The  invariance  of  the  results  from 
shoal  to  shoal  and  day  to  day  is  consistent  with  the  low 
measured  standard  deviations  obtained  for  each  shoal  by  sta¬ 
tionary  averaging.  The  very  strong  roll-off  in  frequency  is 

consistent  with  that  found  just  below  the  resonance  peak  of  a 

38 

system  undergoing  damped  harmonic  oscillation. 

C.  Using  measured  low-frequency  target  strength  to 
infer  swimbladder  properties 

Air-filled  swimbladders  typically  comprise  roughly  5% 
of  fish  body  volume  at  neutral  buoyancy  depth.  To  remain 
neutrally  buoyant  as  hydrostatic  pressure  changes  with 
depth,  fish  need  to  regulate  the  amount  of  gas  in  their  swim- 
bladders  to  maintain  the  5%  ratio.40  Given  this  ratio  and  total 


950  versus  41 5  950  versus  735  1125  versus  950 


FIG.  7.  OAWRS  scattering  strength  level  differences  for  the  indicated  frequency  pairs  as  a  function  of  areal  fish  density  in  dB  re  1  m2  for  data  acquired 
between  22:00  and  22:45  EDT  on  October  3.  The  scattering  strength  difference  at  high  fish  densities  equals  the  target  strength  difference  for  the  frequency  pair 
shown. 
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fish  volume,  swimbladder  volume  can  be  determined  at  any 
depth  as  can  neutral  buoyancy  depth  if  the  relationship  be¬ 
tween  pressure  and  volume  is  known  for  the  swimbladder. 
One  highly  plausible  relationship  is  Boyle’s  law,41'42  where 
the  product  of  pressure  and  volume  remains  constant,  which 
has  been  demonstrated  in  the  laboratory  with  a  single  dead 
herring,43  but  has  not  been  directly  confirmed  in  the  wild 
where  it  is  difficult  to  make  in  situ  measurements  of  the 
physiology  of  free-ranging  fish  at  depth. 

Herring  are  physostomes,  fish  with  open  swimbladders 
connected  to  the  gut  and  colon  44  There  are  three  hypotheses 
by  which  herring  inflate  their  swimbladders:  (1)  by  gulping 
air  at  the  surface,44  47  (2)  by  bacterial  fermentation  in  the 
gut,44'48  and  (3)  by  secretion  of  gas  from  the  blood  stream 
into  the  swimbladder . 49-52  Npttestad5-  found  hypotheses  (1) 
and  (2)  implausible  for  the  Norwegian  spring  spawing  her¬ 
ring  in  his  study  that  remain  in  deep  layers,  are  not  observed 
near  the  surface,  and  are  not  feeding.  This  led  him  to  suggest 
hypothesis  (3).  Nero  et  al.1  conducted  on  site  experiments  on 
the  northern  flank  of  Georges  Bank  with  spawning  herring 
by  adding  weights  to  captured  herring  until  they  sank.  They 
concluded  “that  these  herring  contained  up  to  at  least  a  three 
times  greater  volume  of  gas  than  a  neutrally  buoyant  fish  at 
the  sea  surface,”  and  arrived  at  neutral  buoyancy  depths  as 
great  as  60  m  from  these  and  low-frequency  acoustic  target 
strength  measurements.  They  found  their  results  to  be  con¬ 
sistent  with  measurements  of  Pacific  herring  in  Puget 
sound,48  and  suggested  hypotheses  (2)  and  (3)  as  plausible 
explanations  for  their  observations.  Similarly,  Fange51  sup¬ 
ported  hypothesis  (3)  by  noting  that  “The  herring  lacks  ob¬ 
vious  gas  depository  structures  (rete  mirabile ,  gas  gland),  but 
has  relatively  high  02  values  in  the  swimbladder  (up  to 
21. 5%), 53  and  observations  of  release  of  gas  bubbles  from 
vertically  migrating  herring54  indicate  that  some  gas  secre¬ 
tion  may  occur.” 

The  conditions  of  the  present  experiment  were  not  only 
similar  to  those  of  Nero  et  al.,  but  also  to  those  of 
Npttestad  "  in  that  the  spawning  herring  were  only  observed 
in  deep  layers  and  not  near  the  surface  where  it  is  unsafe  due 
to  predator  attack  as  both  Npttestad  "  and  Makris  et  al. 
noted,  and  a  vast  majority  of  the  herring  captured  in  trawl 
samples  were  observed  to  have  no  large  prey  (copepods)  in 
their  stomachs.  The  latter  point,  however,  needs  to  be  tem¬ 
pered  because  bacterial  content  was  indeterminate  in  the 
samples,  food  resources  are  plentiful  near  the  seafloor  where 
the  herring  shoals  of  the  present  study  were  found,  and  gas 
production  by  bacteria  in  herring  stomachs  can  last  more 
than  90  h  (Ref.  55)  after  ingestion.  It  is  typically  associated 
with  a  delay  due  to  phase  lag  in  bacterial  growth.85  Since  the 
diffusion  rate  of  gas  out  of  the  swimbladders  of  caged  her¬ 
ring  at  fixed  depth  is  also  found  to  be  small,40  corresponding 
to  less  than  a  0.3  dB  decrease  in  target  strength  per  day, 
herring  may  maintain  bacterial  gas  for  long  periods  with 
minimal  feeding  on  large  prey.  Such  feeding  is  known  to 
increase  with  gonadal  development.55  These  facts  suggest 
that  in  addition  to  Npttestad’s52  hypothesis  (3),  hypothesis 
(2)  may  also  remain  highly  plausible  for  our  experiment  as 
suggested  by  Nero  et  al.1  for  a  similar  location  and  season. 

Here  we  compare  the  estimated  low-frequency  target 

114  J.  Acoust.  Soc.  Am.,  Vol.  127,  No.  1,  January  2010 


FIG.  8.  Fork  length  distributions  of  most  frequently  caught  species,  Atlantic 
herring,  Acadian  redfish,  Fladdock,  and  Silver  hake,  from  trawls  deployed 
on  Georges  Bank  (Fig.  1).  The  mean  fork  length  of  herring  is  24.2  cm  with 
a  standard  deviation  6.8%  of  the  24.2  cm  mean.  The  equation  Ltl 
=  1.1  03/.H  +0.0  I  (Ref.  7)  is  used  to  convert  herring’s  fork  length  to  the  total 
length,  where  Ltl  and  Lfl  are  in  cm.  The  mean  fork  length  of  redfish  is 
26.2  cm  with  a  standard  deviation  15%  of  the  26.2  cm  mean.  The  equation 
/.ii  =  1  .033/.h  -0.038  (Ref.  72)  is  used  to  convert  redfish’s  fork  length  to 
total  length.  Silver  hake’s  fork  length  is  the  total  length. 

strength  obtained  from  experimental  data  with  that  derived 
from  Love’s  widely  used  model  for  resonant  scattering  from 
a  fish  swimbladder.6'7'57  This  comparison  enables  estimates 
of  swimbladder  volume  to  be  inferred  for  the  shoaling  her¬ 
ring  observed  in  this  study.  Love  models  the  fish  swimblad¬ 
der  as  an  elongated-spheroidal,  viscous,  heat-conducting 
shell  which  encloses  an  air  cavity  with  surface  tension  at  the 
inner  surface.6’57  This  leads  to  well-understood  damped  reso¬ 
nance  behavior.  While  the  material  in  and  around  the  swim¬ 
bladder  has  more  complex  elastic  composition  and  structural 
constraints  than  that  in  the  Love  model  that  could  lead  to 
more  complex  scattering,  the  Love  model  has  been  success¬ 
fully  tested  in  experimental  settings  where  strong  resonances 
have  been  observed,6'7’57-61  and  probably  provides  an  accu¬ 
rate  description  of  the  first  order  physics  near  resonance. 
Following  Refs.  7  and  61-64,  the  swimbladder  is  modeled  as 
a  resonant,  air-filled  prolate-spheroid  with  a  major-to-minor- 
axis  ratio  of  5-10  and  major-axis  to  fish  length  ratio  of  0.13- 
0.17  based  on  our  trawl  samples  and  CFFS  measurements  for 
herring.  Only  changes  in  minor-axis  are  assumed  to  contrib¬ 
ute  to  swimbladder  volume  change  due  to  physiological 
constraints.7’45’65’66  We  use  the  herring  length  distribution 
[Fig.  8(A)]  measured  from  trawl  samples  and  depth  distribu¬ 
tions  determined  by  CFFS  line-transects.  The  fish  weight 
(W)  to  length  (L)  relationship  is  approximated  by  a  normal 
distribution  with  a  mean  given  by  an  empirically  determined 
length-weight  regression  (gray  line)  from  length-weight 
measurement  of  1219  herring  samples  and  a  standard  devia¬ 
tion  of  ±  20%  of  the  mean  calculated  from  the  length- weight 
data  (dots),  as  shown  in  Fig.  9.  Properties  of  modeled  fish, 
such  as  the  flesh  density,  viscosity,  and  swimbladder  volume 
at  depth  are  tabulated  in  Table  III.  Given  these  constraints, 
swimbladder  volume  or  equivalently  swimbladder  semi- 
minor-axis  is  the  only  unknown  variable  in  the  Love  model 
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FIG.  9.  Atlantic  herring  length-weight  regression  calibration.  The  dots  are 
the  length-weight  data  obtained  from  the  trawl-survey  conducted  by  U.S. 
National  Marine  Fisheries  Service  in  conjunction  with  OAWRS  2006  ex¬ 
periment,  and  the  gray  solid  line  indicates  the  derived  best-fit  length-weight 
regression,  which  can  be  expressed  as  W=aLb,  where  W  is  the  weight  of 
herring  in  kg,  L  is  the  fork  length  of  herring  in  cm,  and  a  and  b  are  empirical 
regression  parameters.  For  this  trawl  dataset,  u  =  3.35X  10-6  and  b— 3.35. 

that  can  lead  to  a  significant  change  in  fish  target  strength  at 
and  below  the  resonance  peak.  Swimbladder  volume  is  as¬ 
sumed  to  vary  with  pressure  according  to  Boyle’s  law,41’42 
from  which  neutral  buoyancy  depth  can  be  uniquely  deter¬ 
mined.  Neutral  buoyancy  depth  is  then  assumed  to  be  a 
Gaussian  random  variable  with  mean  and  standard  deviation 
determined  by  least-squares  fit  between  measured  and  mod¬ 
eled  target  strengths. 

The  best  least  square  fits  between  our  measured  target 
strength  data  and  the  Love  model  appear  in  Figs. 
10(A)-10(E)  for  five  fish  shoals  imaged  by  OAWRS  on  four 
different  days,  and  consistently  show  a  broad  resonance  peak 
with  maximum  at  roughly  1.5 -1.7  kHz  and  swimbladder 
semi-minor-axis  of  3-5  mm.  The  model  to  data  match  is 
excellent,  with  negligible  mean-squared  error,  which  is  sig¬ 
nificant  because  it  would  not  be  possible  to  obtain  a  good 
match  between  the  measured  frequency  dependence  and  the 
Love  model  if  the  overall  level  of  the  measured  target 
strength  data  had  a  significant  scale  factor  error  that  was 
much  larger  than  measured  errors  of  roughly  1  dB.  This  con¬ 
sistency  gives  added  confidence  to  the  veracity  of  both  the 
data  and  model.  As  expected  from  the  roughly  20  dB/octave 
frequency  roll-off  of  the  data,  the  best  fit  of  the  model  is  for 
a  resonance  just  above  the  highest  frequency  data  point 
available  in  the  current  set  of  measurements.  The  model  fits 
of  Fig.  10  can  be  interpreted  with  the  aid  of  Table  VI  which 
shows  the  volume  and  corresponding  swimbladder  minor 
axes  radii  given  the  measured  herring  length  distributions,  as 
well  as  possible  neutral  buoyancy  depths  based  on  Boyle’s 
law.  Neutral  buoyancy  depths  were  found  to  correspond  to 
roughly  half  the  mean  shoal  depth  given  the  measured  her¬ 
ring  length  distribution,  spanning  20-34  cm  with  a  mean  of 
26.7  cm,  and  depth  distribution  in  a  layer  between  120  and 
190  m  from  trawl  and  CFFS  sampling. 

The  localized  measurements  of  Nero  et  al.1  for  herring 
target  strength  data  in  the  1.5 -5.0  kHz  range  show  target 
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strength  levels  consistent  with  those  found  in  our  best-fit 
curve  for  frequencies  above  roughly  2.2  kHz.  This  can  be 
seen  by  noting  their  best-fit  neutral  buoyancy  depth  curve 
(50  m)  falls  within  1-2  dB  of  ours  for  all  examples  above 
roughly  2.2  kHz  for  the  measured  fish  depth,  length,  and 
density  distributions  of  this  study.  The  Love-model  based 
neutral  buoyancy  depths  and  resonances  found  here  are  also 
consistent  with  those  measured  by  OAWRS  in  2003  south  of 
Long  Island,  NY,  in  scattering  from  shoals  that  evidence  sug¬ 
gest  were  also  Atlantic  herring.-'67  Arbitrarily  constraining 
the  neutral  buoyancy  depth  to  be  near  the  surface,  say,  at  4  m 
and  using  Boyle’s  law  for  fish  at  120-190  m  depth  yields  a 
corresponding  3.7  kHz  resonance  and  leads  the  Love  model 
to  a  herring  target  strength  20-30  dB  lower  (Fig.  11)  than 
that  measured  by  the  OAWRS  system.  Such  low  target 
strengths  are  not  only  inconsistent  with  the  target  strength 
measurements  of  Fig.  10  and  Nero  et  al.,  they  also  violate 
conservation  of  energy  since  the  corresponding  transmission 
losses  required  to  match  our  measured  sound  pressure  levels 
returned  from  the  fish  shoals  would  have  to  be  less  than 
those  found  in  even  a  perfectly  reflecting  waveguide  without 
any  medium  attenuation,  i.e.,  the  waveguide  would  have  to 
somehow  add  energy  to  the  source  signal  by  two  to  three 
orders  of  magnitude.  This  can  be  seen  in  Fig.  14  which 
shows  the  measured  transmission  loss  to  be  within  2-3  dB 
of  that  found  in  a  perfectly  reflecting  waveguide  without 
attenuation. 

A  wide  distribution  of  swimbladder  volumes  and  corre¬ 
sponding  neutral  buoyancy  depths  within  any  shoal  is  likely 
and  could  potentially  unify  the  various  existing  data  sets  by 
superposition,  with  larger  swimbladder  volumes  dominating 
at  the  lower  end  of  the  possible  resonance  spectrums  and 
smaller  volumes  dominating  the  higher  end.  None  of  the  sys¬ 
tems  used  in  the  field  so  far,  however,  could  test  this  since  it 
would  require  acquisition  of  simultaneous  data  both  well  be¬ 
low  and  well  above  all  contributing  resonant  frequencies. 
Since  low-frequency  target  strength  measurements  near  reso¬ 
nance  are  far  more  sensitive  to  small  changes  in  swimbladder 
minor-axis  or  volume  than  at  much  higher  CFFS  frequencies, 
such  lower  frequency  measurements  have  the  potential  to 
resolve  the  in  situ  swimbladder  volumes  of  wild  herring  at 
depth  with  much  greater  accuracy.  The  measurements  of  tar¬ 
get  strength  for  herring  as  a  function  of  depth  reported  by 
Ona  at  38  kHz,  for  example,  show  standard  deviations  of 
8-10  dB  and  minimal  3-5  dB  reduction  from  the  surface  to 
300  m,66  which  falls  well  within  these  standard  deviations. 
These  standard  deviations  significantly  exceed  the  5.6  dB  ex¬ 
pected  purely  from  stationary  Gaussian  field  fluctuations  by 
the  Central  limit  theorem,-6  shown  to  apply  to  scattering 
from  fish  of  random  orientation  by  Dahl  and  Mathisen.68 
This  may  suggest  a  bias  in  fish  orientations  or  a  wide  spread 
in  swimbladder  volumes  in  the  Ona  data,  which  could  eas¬ 
ily  include  the  span  of  volumes  measured  here  and  in  Ref.  7. 
Fassler  et  a/.43  suggested  that  projected  dorsal  area  variations 
exhibit  greater  depth  variation  than  the  mean  target  strength 
measurements  of  Ona.'  We  find  the  Ona'  mean  trend  in 
depth  to  be  consistent  with  that  expected  from  Boyle’s  law 
since  at  their  frequencies  ka,  the  product  of  acoustic  wave- 
number  to  the  semi-minor-axis  of  the  swimbladder,  is  typi- 
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FIG.  10.  Experimentally  determined  low-frequency  target  strength  corresponding  to  the  average  scattering  cross-sections  of  shoaling  herring  observed  from 
OAWRS  imagery  acquired  from  five  shoals  on  4  days  (Table  VI)  at  415,  735,  950,  and  1125  (circles)  with  standard  deviations  (error  bars).  Comparison  with 
Love-model  mean  target  strength  for  shoaling  herring,  with  physical  parameters  tabulated  in  Table  III,  of  different  swimbladder  semi-minor  axes  over  the 
shoals’  depth  distributions  (lines).  The  best  least-squares  fits  shown  are  obtained  only  using  target  strength  estimates  of  each  shoal  (Table  IV)  with  standard 
deviations  less  than  3  dB.  Arrows  indicate  the  target  strength  uncertainties  due  to  potential  masking  from  background  scattering  (Sec.  HIE)  for  given 
frequencies.  The  best-fit  means  and  standard  deviations  of  inferred  swimbladder  volume,  swimbladder  semi-minor  axes,  and  corresponding  neutral  buoyancy 
depths  of  each  shoal  were  tabulated  in  Table  VI.  (A)  Shoaling  herring,  distributed  between  120  and  185  m  (Fig.  4),  imaged  with  the  OAWRS  system  from 
23:30  EDT  October  2  to  01:00  EDT  on  October  3.  (B)  Shoaling  herring,  distributed  between  135  and  175  m,  imaged  with  the  OAWRS  system  from  18:55 
to  19:50  EDT  on  October  3.  (C)  Shoaling  herring,  distributed  between  120  and  175  m,  imaged  with  the  OAWRS  system  from  22:00  to  22:45  EDT  on  October 
3.  (D)  Shoaling  herring,  distributed  between  120  and  185  m,  imaged  with  the  OAWRS  system  from  06:10  to  06:50  EDT  on  September  27  (Fig.  13).  (E) 
Shoaling  herring,  distributed  between  150  and  180  m,  imaged  with  the  OAWRS  system  from  07:25  to  07:50  EDT  on  September  29. 


cally  less  than  unity  and  target  strength  is  no  longer  propor¬ 
tional  to  projected  area  as  in  the  large  ka  limit. 

Evidence  from  the  extensive  OAWRS,  CFFS,  and  trawl 
surveys  conducted  during  OAWRS  2006  offers  no  plausible 
alternative  to  herring  as  the  primary  constituent  and  source 
of  scattering  in  the  shoals  imaged  by  OAWRS.  Consideration 
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of  the  effect  of  two  other  species,  however,  also  present  in 
many  trawls  but  in  far  fewer  numbers  (Table  VII),  still  pro¬ 
vides  useful  perspective.  These  are  silver  hake  and  Acadian 
redfish.  While  silver  hake  were  found  at  shallower  depths 
than  herring,  their  much  shorter  lengths  (Fig.  8)  lead  to  in¬ 
dividual  target  strengths  more  than  20  dB  lower  than  those 
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TABLE  VI.  Experimentally  inferred  means  and  standard  deviations  of  swimbladder  volume,  vzh  and  <tv^, 
semi- minor-axis,  dz  and  <ra  ,  over  the  depth  distributions  of  the  shoals,  and  corresponding  means  and  standard 
deviations  of  neutral  buoyancy  depth,  dnh  and  where  neutral  buoyancy  depth  is  restricted  to  water-column 
depths  of  0-200  m  in  the  least  squares  fit.  All  three  parameters  are  assumed  to  be  Gaussian  random  variables 
completely  characterized  by  their  respective  means  and  standard  deviations. 
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“Time  periods  during  which  both  OAWRS  and  CFFS  systems  simultaneously  co-registered  contiguous  shoal 
segments.  The  time  is  in  EDT. 


measured  at  the  location  of  the  herring  shoals  at  OAWRS 
frequencies  (Fig.  11).  Given  that  (1)  both  CFFS  and  trawl 
samples  found  silver  hake  in  considerably  lower  areal  densi¬ 
ties  than  herring,  (2)  their  individual  target  strengths  are 
roughly  20  dB  lower,  and  (3)  OAWRS  transmission  loss  was 
much  greater  in  the  shallower  layers  where  silver  hake  re¬ 
sided  as  part  of  the  experimental  design,  contributions  from 
silver  hake  can  be  ruled  out  as  a  plausible  explanation  for 
shoals  imaged  by  OAWRS.  As  shown  in  Fig.  11(A),  the  pres¬ 
ence  of  redfish  in  the  maximum  percentages  determined  by 
trawl  has  a  negligible  effect  on  the  best-fit  target  strength 
compared  to  that  of  modeling  herring  alone  [Figs. 
10(A)-10(E)].  This  should  be  expected  because  the  low- 
frequency  dependence  of  background  returns  when  no  shoals 
are  present  is  effectively  negligible  compared  to  the  depen¬ 
dence  when  shoals  are  present  (Fig.  5)  and  simultaneous 
trawls  showed  the  shoals  to  be  overwhelmingly  comprised  of 
herring  but  the  background  to  yield  negligible  amounts  of 


herring  (Table  IV,  Fig.  12).  So,  even  though  redfish  are 
physoclists  with  trawl  sample  lengths  (Fig.  8)  typically 
greater  than  those  of  the  herring,  the  measured  areal  densities 
of  redfish  are  too  low  to  have  a  significant  impact  on  the 
average  scattering  cross-section  measured  in  the  observed 
shoals.  Including  greater  percentages  of  redfish  in  the  mod¬ 
eling  or  including  measured  percentages  and  shallower  neu¬ 
tral  buoyancy  depths  for  herring  leads  to  far  worse  fits  and 
high  mean-square  error  when  matched  with  the  data,  as 
shown  in  Fig.  11(B). 

D.  Space,  time,  and  frequency  dependencies 

From  another  perspective,  since  the  spatio-temporal 
population  distributions  of  the  large  shoals  versus  back¬ 
ground  levels  are  consistent  among  OAWRS,  CFFS,  and 
trawls,  as  shown  here  and  in  Ref.  3,  the  spatio-temporal  dis¬ 
tributions  of  silver  hake,  redfish,  or  any  other  contaminant 
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FIG.  1 1 .  (A)  An  example  of  Love-model  target  strength  corresponding  to  the  average  scattering  cross-section  of  an  individual  for  mixed  species  content  and 
swimbladder  semi-minor  axes  shown.  Comparison  of  Love-model  target  strength  with  experimentally  determined  mean  target  strength  estimates  of  shoaling 
herring,  distributed  between  135  and  175  m,  imaged  with  the  OAWRS  system  from  18:55  to  19:50  EDT  on  October  3.  Presence  of  redfish  in  trawl-determined 
percentage  (dashed  gray  line)  has  negligible  effect  on  best-fit  target  strength  compared  to  herring  alone  (solid  black  line),  while  including  an  unrealistic 
percentage  of  redfish  yields  far  worse  fits  (dashed  black  line).  The  silver  hake,  found  at  shallower  water  depth  (<75  m)  in  trawls,  have  resonance  peak  above 
3  kHz  (dash-dot  curve)  making  their  contribution  negligible.  Herring  target  strength  with  a  resonance  frequency  at  3.7  kHz  (dashed  gray  line)  based  on  Love’s 
model  using  length  and  depth  distribution  obtained  from  CFFS  and  trawl  measurements  is  found  to  be  neutrally  buoyant  at  4  m,  and  is  20-30  dB  lower  than 
those  measured  by  the  OAWRS  system.  (B)  Same  as  (A)  but  for  shallower  herring  neutral  buoyancy  depths.  Arrows  indicate  potential  target  strength 
uncertainties  for  given  frequencies. 
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TABLE  VII.  Concurrent  high-speed  rope  trawl  deployed  by  NOAA  FRV  Delaware  II  within  or  in  the  vicinity 
of  large  herring  shoals  imaged  by  OAWRS  during  the  OAWRS  2006  experiment  in  the  Gulf  of  Maine  and 
Georges  Bank  at  shoal  depth.  The  number  of  most  frequently  caught  species  in  each  trawl  deployment,  includ¬ 
ing  Atlantic  herring,  Acadian  redfish,  Silver  hake,  and  Haddock  are  tabulated.  Trawls  134,  137,  and  139  were 
made  with  simultaneous  OAWRS  imagery,  and  trawls  137  and  139  were  made  directly  through  shoals  as  shown 
in  Figs.  12(B)  and  12(C).  In  contrast,  trawl  134  was  made  in  a  region  with  no  shoal  formed  (Fig.  12(A)),  but 
one  shoal  would  form  in  the  vicinity  4  h  later.  Trawl  105  was  made  through  shoals  imaged  by  OAWRS  95  min 
before,  and  trawl  106  was  made  through  shoals  imaged  by  OAWRS  20  min  later. 


Deployment 

Atlantic  herring 

Acadian  redfish 

Silver  hake 

Haddock 

Total  catch 

105 

8030  (99.98%) 

0 

0 

0 

8032 

106 

634  (96.79%) 

0 

14  (2.14%) 

2  (0.31%) 

655 

134 

3  (11.54%) 

0 

14  (53.85%) 

0 

26 

137 

333  (76.03 ) 

0 

94  (21.46%) 

0 

438 

139 

796  (74.74%) 

9  (0.85%) 

208  (19.53%) 

23  (2.16%) 

1065 

Total  catch  of  experiment 

9796  (96.54%) 

9  (0.07%) 

330  (2.67%) 

25  (0.2%) 

10216 

species  would  have  to  consistently  follow  those  of  herring  if 
they  were  a  major  contributor  to  OAWRS  returns,  which  is 
both  implausible  and  contrary  to  trawl  and  CFFS  data  as  well 
as  the  frequency  dependence  of  OAWRS  returns  in  shoals 
versus  background  (Fig.  5).  Also,  neither  our  trawl  nor  CFFS 
data  show  any  evidence  of  shallow  fish  layers  that  could 
account  for  the  prominent  OAWRS  returns  that  co-registered 
in  space  and  time  with  the  deep  shoals  measured  by  simul¬ 
taneous  trawl  and  CFFS  sampling.  Indeed,  considering  ex¬ 
perimental  causality,  it  was  consistently  necessary  for 
OAWRS  to  first  find  and  then  direct  trawl  and  CFFS  vessels 
to  the  locations  of  these  deep  shoals  because  they  are  so 
difficult  to  find  with  conventional  methods  given  the  fact 
they  occupy  areas  many  orders  of  magnitude  smaller  than  the 
wide  areas  over  which  they  may  roam. 

V.  CONCLUSION 

The  low-frequency  target  strength  of  Atlantic  herring 
(' Clupea  harengus )  is  estimated  from  experimental  data  ac¬ 
quired  from  shoaling  herring  in  the  Gulf  of  Maine  during  the 


Autumn  2006  spawning  season  in  the  300-1200  Hz  range 
using  simultaneous  ocean  acoustic  waveguide  remote  sens¬ 
ing,  conventional  fish  finding  sonar,  and  trawl  surveys.  The 
target  strength  expected  of  an  individual  is  found  to  have  a 
strong  nonlinear  dependence  on  frequency  consistent  with 
resonant  scattering  from  an  air-filled  swimbladder  given 
measured  fish  length,  depth  distributions,  and  experimentally 
inferred  swimbladder  volumes  based  on  Love’s  model, 
which  indicate  the  herring  remain  negatively  buoyant  in  lay¬ 
ers  near  the  seafloor  for  extended  periods.  The  OAWRS  sys¬ 
tem  used  in  this  study  employed  an  instantaneous  imaging 
diameter  of  100  km  with  regular  minute-to-minute  updates 
enabling  unaliased  monitoring  of  fish  populations  over  eco¬ 
system  scales.  This  included  detection  and  imaging  of  shoals 
of  Atlantic  herring  containing  hundreds  of  millions  of  indi¬ 
viduals,  as  confirmed  by  concurrent  trawl  and  CFFS  surveys 
that  were  directed  to  the  shoals’  locations  by  OAWRS.  High 
spatial-temporal  coregistration  was  found  between  herring 
shoals  imaged  by  OAWRS  and  concurrent  CFFS  line- 
transects. 


Eastings  to 
OAWRS  source  (km) 

FIG.  1 2.  (Color  online)  [(A)-(C)]  Locations  of  trawls  over  simultaneous  OAWRS  images.  Trawls  137  and  139  were  made  directly  through  shoals  as  shown 
in  (B)  and  (C).  In  contrast,  trawl  134  was  made  in  a  region  with  no  shoal  (A),  but  one  shoal  would  form  in  the  vicinity  4  h  later.  The  OAWRS  source  locations 
are  the  coordinate  origin  in  all  OAWRS  images.  On  October  3,  the  OAWRS  source  ship  was  moored  at  42.2089N,  67.6892W.  The  trawls  in  (A)-(C)  were 
deployed  and  towed  along  the  solid  lines  starting  at  a  and  ending  at  Cl.  The  dashed  lines  indicate  the  contours  of  100,  150,  180,  and  200  m  water  depth. 
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TABLE  VIII.  Concurrent  high-speed  rope  trawl  deployed  by  NOAA  ship  FRV  Delaware  II  within  or  in  the 
vicinity  of  large  herring  shoals  imaged  by  OAWRS  system  during  the  OAWRS  2006  experiment  in  the  Gulf  of 
Maine  and  Georges  Bank  at  shoal  depth.  The  dates,  times  (Eastern  Daylight  Time),  deploy  depths,  and  geo¬ 
graphic  locations  of  the  trawls  are  tabulated. 


Deployment 

Date 

Time 

Deploy  depth 
(m) 

Begin  Lat 

Begin  Lon 

End  Lat 

End  Lon 

105 

09/26 

14:21:40-15:16:32 

130-160 

41°  53.00N 

68°  06.62W 

41°  51.95N 

68°  10.45W 

106 

09/26 

19:23:24-19:57:20 

130-160 

41°  55.49N 

68°  03.80W 

41°  55.05N 

68°  06.33W 

134 

10/03 

13:21:12-13:54:01 

140-170 

42°  01.80N 

67°  53.39W 

42°  03.06N 

67°  55.24W 

137 

10/03 

18:47:09-19:27:04 

140-170 

42°  06.36N 

67°  41.02W 

42°  04.76N 

67°  43.27W 

139 

10/03 

21:16:35-21:52:46 

140-170 

42°  08.54N 

67°  33.83W 

42°  07.00N 

67°  35.63W 
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APPENDIX  A:  CONCURRENT  NATIONAL  MARINE 
FISHERIES  SERVICE  ANNUAL  ATLANTIC  HERRING 
ACOUSTIC  SURVEY  IN  THE  GULF  OF  MAINE 

In  conjunction  with  OAWRS  2006  experiment,  U.S.  Na¬ 
tional  Marine  Fisheries  Service  conducted  annual  Atlantic 
herring  survey  of  the  Gulf  of  Maine  and  Georges  Bank.  Five 
trawls,  105,  106,  134,  137,  and  139,  targeted  at  130-170  m 
depth  (Table  VIII)  were  deployed  within  or  in  the  vicinity  of 
the  large  fish  shoals  imaged  by  OAWRS  system,  which  en¬ 
abled  onsite  species  identification  and  biological  measure¬ 
ments.  The  locations  of  these  trawls  are  shown  in  Fig.  1.  The 
trawls  134,  137,  and  139  were  made  with  simultaneous 
OAWRS  imagery,  and  trawls  137  and  139  were  made  di¬ 
rectly  through  shoals,  as  shown  in  Figs.  12(B)  and  12(C).  As 
tabulated  in  Table  VII,  Atlantic  herring  is  found  consistently 
to  be  the  overwhelmingly  predominant  species  in  trawls  105, 
106,  137,  and  139.  Trawl  134,  shown  in  Fig.  12(A),  was 
deployed  in  a  region  with  no  shoal  formed,  but  one  would 
form  in  the  vicinity  4  h  later,  and  like  both  OAWRS  and 
CFFS  does  not  show  high  concentrations  of  herring  before 
the  shoal  formed.  Trawl  105  was  made  through  shoals  im¬ 
aged  by  OAWRS  95  min  before,  and  trawl  106  was  made 
through  shoals  imaged  by  OAWRS  20  min  later.  In  addition, 
a  small  percentage  of  juvenile  silver  hake,  Acadian  redfish, 
and  haddock  are  also  present  in  the  trawls.  Most  of  the 
juvenile  silver  hake  were  caught  at  shallower  water  depth 
(<75  m)  as  the  deeper  trawls  were  deployed.  Histograms  of 
the  measured  length  distributions  of  the  most  frequently 
caught  species  are  plotted  in  Fig.  8.  The  conversion  of  her¬ 
ring  length  to  weight  by  regression  analysis  is  given  by  W 
=  0.003  35L3'35,  based  on  length-weight  measurements  of 
1219  herring  caught  in  trawls  throughout  the  course  of  the 
survey. 
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APPENDIX  B:  ESTIMATING  LOW-FREQUENCY 
TARGET  STRENGTH  OF  SHOALING  HERRING 
WITHOUT  CALIBRATED  TARGET  DEPLOYMENT 

Here  we  present  the  estimated  low-frequency  target 
strength  of  a  shoaling  herring  based  on  simultaneous 
OAWRS  and  CFFS  sampling  line-transects  of  large  shoals 
on  September  27  and  29,  2006  when  no  calibrated  targets 
were  available  to  aid  in  charting.  Only  minor  adjustments  of 
less  than  1.5°  from  the  array  heading  sensors  were  needed  to 
ensure  optimal  coregistration  between  OAWRS  and  CFFS 
locations  of  shoals,  which  were  typically  the  only  landmark 
available  in  the  absence  of  the  calibrated  targets.  The  target 
strength  estimates  shown  here  are  consistent  with  those  pre¬ 
sented  for  the  other  days  when  calibrated  targets  were  avail¬ 
able  in  the  main  text. 

During  the  early  morning  hours  of  September  27,  be¬ 
tween  06:10  and  06:50  EDT,  both  OAWRS  and  CFFS  co¬ 
registered  a  massive  herring  shoal  10-15  km  southeast  of 
the  OAWRS  source,  spanning  more  than  2  km  in  range  and 
5  km  in  azimuth,  as  shown  in  Figs.  13(A)— 13(C).  A  dense 
layer  of  herring,  shown  in  Fig.  13(D),  was  found  consistently 
spanning  120-170  m  in  the  water  column  by  simultaneous 
CFFS  line-transects.  The  resulting  target  strength  estimates 
at  950  Hz  obtained  by  means  of  the  approach  described  in 
Sec.  Ill  D  are  plotted  in  Fig.  13(G).  The  target  strength  esti¬ 
mates  at  415,  735,  and  1125  Hz  are  tabulated  in  Table  IV 
using  the  approach  described  in  Sec.  Ill  E,  where  target 
strength  at  950  Hz  is  used  as  the  base  frequency.  No  conclu¬ 
sion  can  be  made  for  the  target  strength  difference  between 
950  and  415  Hz,  because  fish  density  is  not  high  enough  to 
make  scattering  from  herring  at  415  Hz  distinguishable  from 
background  levels.  The  same  approaches  are  also  applied  to 
estimate  the  target  strength  expected  for  an  individual  shoal¬ 
ing  herring  imaged  on  September  29  between  07:25  and 
07:50  EDT.  The  results  are  provided  in  Table  IV. 


APPENDIX  C:  STATISTICAL  ANALYSIS  FOR  LOW- 
FREQUENCY  TARGET  STRENGTH  OF 
HERRING 

Here  we  describe  the  statistical  approach  used  to  esti¬ 
mate  the  mean  and  standard  deviation  of  the  low-frequency 
target  strength  of  fish  populations  imaged  concurrently  by 
OAWRS  and  CFFS.  We  confine  our  analysis  to  regions 
where  both  systems  image  continuous  segments  of  statisti- 
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FIG.  13.  (Color  online)  Herring  target  strength  at  950  Hz  estimated  by  matching  the  areal  density  of  OAWRS  and  CFFS  data  acquired  during  the  early 
morning  hours  of  September  27.  Similar  to  Fig.  4  but  for  a  contiguous  shoal  40-min  segment  starting  at  06:10  EDT  on  September  27. 


cally  stationary  fish  populations  that  span  many  resolution 
cells  in  both  the  CFFS  and  OAWRS  systems. 

Employing  Eq.  (5),  let  TS,-  be  independent  estimates  of 
low-frequency  target  strength  obtained  for  each  resolution 
cell  after  sufficient  spatial  and  temporal  averaging  in  anti-log 
to  eliminate  the  bias  when  converted  to  log  units  and  to 
insure  that  the  TS,  are  Gaussian  random  variables."6  Assum¬ 
ing  stationarity,  let  their  means  be  TS  and  standard  devia¬ 
tions  be  crxs.  Then,  given  N  independent  measurements  of 
the  target  strength  along  the  CFFS  transect,  the  resulting  lin¬ 
ear  estimator  for  mean  target  strength  is69 


°Ts  - 


(N-  1) 


2  (TS,-  -  TS)2. 


(C2) 


The  variance  of  the  mean  target  strength  estimate  is 

Var(TS)  =  — . 

N 


(C3) 


Applying  Eq.  (C3),  the  standard  deviation  in  the  mean  target 
strength  estimate  for  the  data  in  Fig.  4  is  roughly  on  the  order 


of  0.7-0. 8  dB. 


ts  =  -Sts,-  (Cl) 

N ti 

and  the  corresponding  estimator  for  the  standard  deviation  of 
the  target  strength  distribution  is  obtained  by69 
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APPENDIX  D:  QUANTIFYING  ERROR  IN  ABUNDANCE 
ESTIMATES 

Here  we  quantify  the  error  in  the  abundance  estimates 
from  OAWRS  data  presented  in  this  paper.  From  Eqs.  (4) 
and  (5),  error  in  «a,oawrs  at  any  pixel  is  caused  by  a  fixed 
OAWRS  target  strength  estimate,  which  is  the  same  for  all 
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Calibrated  mean  TL  with  instantaneous  MF  data,  (925-975  Hz)  Calibrated  mean  TL  with  6-sample  averaged  MF  data,  (925-975  Hz) 


Range  (km)  Range  (km) 

FIG.  14.  (A)  Experimentally  determined  mean  and  standard  deviation  (4.78  dB)  of  677  measured  instantaneous  broadband  one-way  TL  data  after  matched 
filter  with  950  Hz  center  frequency.  Plotted  as  a  function  of  range,  with  modeled  one-way  TL  overlain.  Transmission  data  acquired  by  a  single  desensitized 
hydrophone  on  OAWRS  receiver  array  on  October  1-3  2006.  Modeled  TL  computed  by  Monte-Carlo  simulation  with  parabolic  equation  inputing  measured 
oceanography  and  bathymetry  of  Georges  Bank  environment.  Modeled  TL  (dashed  gray  line)  in  a  perfectly  reflective  waveguide  with  no  water-column 
attenuation  by  normal-mode  model.  (B)  After  six-sample  averaging,  the  rms  error  is  reduced  to  1.7  dB,  consistent  with  theory  for  stationary  averaging  of 
intensities  of  circular  complex  Gaussian  random  fields  (Refs.  16  and  26). 


pixels,  and  the  OAWRS  scattering  strength,  which  varies 
from  pixel  to  pixel.  Summing  large  numbers  of  pixels  re¬ 
duces  the  percentage  error  from  scattering  strength  fluctua¬ 
tions  in  the  population  density  estimate  to  negligible  values 
by  the  law  of  large  numbers,  leaving  the  error  to  be  domi¬ 
nated  by  that  in  the  OAWRS  target  strength  estimate,  which 
for  the  950  Hz  imagery  is  consistently  between  0.7  and 
0.8  dB  throughout  the  experiment,  corresponding  to  an  error 
of  17-20%  in  total  population. 


APPENDIX  E:  CALIBRATION  OF  BROADBAND  MEAN 
TL  WITH  ONE-WAY  PROPAGATED  EXPERIMENTAL 
DATA  IN  THE  NORTHERN  FLANK  OF 
GEORGES  BANK 

In  wide-area  sonar  applications,  two-way  TL  must  be 
efficiently  and  accurately  estimated  in  order  to  invert  for 
scattering  strength  or  target  strength  of  scatterers.  ’  Here, 
RAM  is  used  to  estimate  two-way  TL  expected  over  wide 
areas  with  measured  bathymetry  and  oceanography  of  the 
northern  flank  of  Georges  Bank  environment.  For  each 
OAWRS  image,  TL  was  computed  along  radials  separated  by 
roughly  1.5°,  half  the  receiver  array  broadside  resolution,  for 
forward  transmission  from  the  source,  and  again  for  return 
transmission  to  the  receiver  array,  by  multiple  Monte-Carlo 
realizations  per  radial.  Each  Monte-Carlo  realization  em¬ 
ployed  a  different  measured  sound  speed  profile  every  500  m 
to  incorporate  the  effects  of  the  fluctuating  ocean  waveguide. 

In  the  OAWRS  2006  experiment,  a  very  complete  set  of 
more  than  12  000  transmission  loss  measurements  spanning 
ranges  from  hundreds  of  meters  to  tens  of  kilometers,  and 
more  than  roughly  200  sound  speed  profile  measurements 
were  made  over  the  survey  area.  These  were  used  to  calibrate 
the  parabolic  equation  model  used  to  estimate  wide-area 
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transmission  loss,  population  density,  and  target  strength. 
Measured  sound  speed  was  found  to  follow  a  roughly  Gauss¬ 
ian  distribution  about  the  mean  at  each  depth  and  to  be  rela¬ 
tively  uniform  in  time  and  space  over  the  entire  OAWRS 
imaging  area  (Fig.  3)  with  no  evidence  of  horizontal  features 
besides  the  expected  short  term  fluctuations  from  internal 
wave  activity. 

The  best  fit  between  measured  and  modeled  transmis¬ 
sion  loss  was  obtained  for  expected  sandy  bottom  conditions 
of  sound  speed  of  1.7  km  s-1,  density  of  1.9  g  cm-3  and  at¬ 
tenuation  of  0.8  dB  X-1,  and  in-water-column  attenuation  of 
6X  10-5  dB  X-1,  as  shown  in  Fig.  14  using  the  maximum 
likelihood  method.15  Figure  14  shows  the  standard  deviation 
of  one-way  TL  for  instantaneous  broadband  measurements  is 
roughly  4.8  dB,  and  after  six-sample  averaging  is  on  the  or¬ 
der  of  1.7  dB,  following  previous  experimental  results  and 
theory.16,26  We  find  the  mean  measured  TL  to  be  within 
2-3  dB  of  that  of  a  perfectly  reflecting  lossless  waveguide 
of  the  same  depth  (dashed  gray  line  in  Fig.  14).  The  scattered 
field  from  large  fish  shoals  after  two-way  propagation  in  a 
waveguide  also  obeys  circular  complex  Gaussian  field  statis¬ 
tics  by  the  central  limit  theorem  and  has  experimentally  mea¬ 
sured  standard  deviation  of  roughly  1.5  dB,  after  our  stan¬ 
dard  six-sample  averaging,  for  the  current  experiment4  which 
also  follows  that  expected  from  theory  and  numerous  past 
experiments  where  internal  waves,  ocean  turbulence,  eddies, 
and  variations  in  scatterer  orientation  have  been  shown  to 
cause  short  term  intensity  fluctuations  with  a  standard  devia¬ 
tion  that  can  be  reduced  to  a  small  fraction  of  the  mean  after 

.  .•  •  20,21,24,26,67,70 

stationary  averaging. 
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